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1. Introduction 

The simplest vortex dynamics in a superfluid corresponds to the two-dimensional motion 
of a rectilinear vortex filament [Tj. In fact, in an infinite superfluid a vortex "charged" 
with one quantum of counterclockwise circulation will move like an electron in a uniform 
magnetic field, ie performing a circular cyclotron motion ruled by: 

m v r = p s hz x f . (1) 

Here m v denotes the vortex effective mass per unit length, p s the number density of 
the background superfluid at rest, h the Planck's constant and r = (x, y) the two- 
dimensional coordinate of the vortex core. We note that the Magnus force in the right- 
hand side of (0) is formally equivalent to the Lorentz force on a negative point charge in 
a uniform magnetic field parallel to the z axis. Actually, this electromagnetic analogy is 
only a part of a whole mapping by which a 2-D homogeneous superfluid can be mapped 
onto a (2+l)-D electrodynamic system, with vortices and phonons playing the role of 
charges and photons, respectively j2]. For instance, any accelerated motion of a vortex 
would result in the radiation of sound waves in the superfluid, ie the emission of phonons, 
a process which is entirely analogous to the photon radiation mechanism stemming from 
an accelerated charge in electrodynamics. In practice, however, this simple picture only 
would apply to a superfluid formed by 4 He atoms, viz the boson isotope of helium, at 
low temperatures (T <0.4 K). Ordinary helium, on the other hand, contains a small 
amount of impurity, fermion 3 He atoms, which however produces a viscous drag force 
on a moving vortex at the lowest temperatures. In fact, below 0.4 K the scattering of 
thermal phonons by the vortex has negligible effects compared to the drag force due 
to 3 He scattering In addition, at higher temperatures (T >0.5 K), most of the 
elementary excitations of the superfluid 4 He that collide with the vortex are rotons, ie 
quasiparticles having their momentum around the minimum of the dispersion curve. In 
fact, (see figure 1) only the elementary excitations with momentum below 0.5 A -1 can 
be regarded as phonons arising from a linear dispersion relation uj = c s k, whereas the 
rest of the dispersion curve do not yield thermal elementary excitations, except for a 
small interval around the minimum (~ 1.9 A -1 ). Quasiparticles with momentum at 
the right (left) of this minimum have their group velocity parallel (antiparallel) to their 
momentum and are called R + {R~) rotons. Above 0.5 K the source of the drag force 
on a vortex in ordinary helium is roton scattering, being the effect of 3 He collisions 
practically negligible. Whatever the source, however, such a drag force can be written 
as two additional terms in the right-hand side of equation ((TJ), viz 

m v r = (p s h — D')z xr-Dr, (2) 

where D' and D denote transversal and longitudinal friction coefficients, respectively 
[U Uj- It is interesting to compare in figure 1 the dispersion curve for elementary 
excitations in superfluid 4 He, with the energy spectrum of 3 He atoms in ordinary helium. 
Such solvated atoms behave like heavier free particles {et = h 2 k 2 /2m*) with an effective 
mass m* which exceeds two times the mass of a bare 3 He atom. As a final remark about 




Figure 1. Dispersion curve for elementary excitations in superfiuid 4 He (full curve) 
and energy spectrum of solvated 3 He atoms in ordinary helium (broken curve). 



figure 1 we note that the terminations of both curves are due to unstabilities caused 
by roton creation processes. That is, any elementary excitation exceeding two times 
the roton energy should be unstable against decay into two rotons, whereas 3 He atoms 
exceeding the roton energy should decay into a low energy atom plus a roton. 

Here it is important also to take into account another consequence of the above 
scattering processes, apart from the friction itself, that is the thermal excitation of 
vortex waves |H] . These are helical waves in which each vortex line element executes 
a circular motion about the undisturbed line (z axis). The radius of such a circle 
is assumed to be much smaller than the wavelength, so the above elements will keep 
almost parallel to the z axis, fulfilling an equation of motion like (£[}. In fact, such an 
equation has to be generalized to the situation where there exists an external superfiuid 
flow of velocity v s : 

m v r = p s hz x (f - v s ). (3) 
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Such an 'external' superflow corresponds in our case to the local self-induced velocity 
generated by the vortex line curvature PIE], v s = —Vi6, which points in a direction 
opposite to the one of the superfluid velocity field generated by the undisturbed vortex 
line. This self-induced velocity, being proportional to the line displacement |r| from the 
z axis, can be written as Vi = ou-\r\, where 

frk 2 

u.{k) ~ -_[ln(|fc|a) + 0.116] (4) 
2m4 

corresponds to the well-known dispersion relation for Kelvin waves of long wavelength 
A, \k\a <C 1, being \k\ = 2tt/X, a ~ 1 A= vortex core parameter and m 4 = mass of 
a 4 He atom. Notice that the wave vector k, which points along the z axis, can take 
positive or negative values depending on the two ways of generating the vortex helix. 
Then, equation (jSJ) can be rewritten 

r = Qzxr — Q u;_r, (5) 

where 

Q = p s h/m v (6) 

corresponds to the cyclotron frequency arising from the equation of motion (JJJ. Actually, 
it is easy to check that the equation (jSJ) allows both possible directions for circular 
motion, since a replacement |r| = const, r = u |r| in (j2J) leads to a quadratic equation 
in the angular frequency uj with solutions: 

n 



w (±) = 




(7) 



That is, in the limit u^/Q C 1 we have either the counterclockwise cyclotron motion 
of frequency ~ fi, or the usual clockwise polarization of Kelvin waves cu^ ~ — u_ 
(cf (jlj)), for vortices of counterclockwise circulation and negligible mass (Q — > oo). 
Actually, there is no experimental data about vortex trajectories, so the value of the 
vortex mass m v and hence of Q, can only be extracted from theoretical considerations. 
If m v is calculated from a classical hydrodynamical model, Q should be about 3 ps~ x [T], 
whereas more recent theories, for which m v should be logarithmically divergent with the 
system size, lead to JTs of order 0.1 — 0.01 ps _1 , for typical experimental conditions [HE]. 
All these figures are consistent with the approximation Q ^> cj_ for long wavelengths 
| A; | a <C 1, where the dependence on k of the frequency u/ + ) can be neglected [I]. 

Finally, the vortex equation of motion is obtained by adding both effects, friction 
(0) and oscillations (0) together: 



n 



D' \ A . D . 

— - \ z x r -r — uj-Y 

Pshl p s h 



(8) 



Note that in addition to the time dependence of r, one should take into account a 
parametric dependence r(z) following the helix curvature. 

To this point we have regarded the vortex coordinates as classical time-dependent 
variables, but it is important to observe in this respect that the above theoretical 
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estimates of the cyclotron frequency yield values of htt/keT greater than ~0.1 for 
T < 1 K. This seems to indicate that a classical treatment cannot be wholly satisfactory. 
Moreover, even in the case of purely low-frequency Kelvin waves, the need of a quantum 
mechanical analysis was early pointed out by Fetter [§] . Such a theory was in fact used 
to study phonon scattering by a vortex [TH) . and it is our purpose to present in this 
article a more general treatment at which the vortex mass, and hence the cyclotron 
frequency, are included in the theory and assigned finite values. Our starting point will 
be a vortex Hamiltonian from which the equation of motion ((HJ) derives. Then, after 
quantization of the vortex variables, we will show that such a Hamiltonian consists of 
independent harmonic oscillator modes of frequencies Q and cj_, which interact with 
the heat bath represented by the ordinary helium at a finite temperature. Such an 
interaction is modelled through a generic momentum-conserving scattering Hamiltonian, 
which is used to study the dissipative dynamics of the cyclotron modes. Thus, we shall 
show that the friction turns out to be essentially unaffected by such oscillations, allowing 
us to extend our previous conclusions on the memory effects on straight vortex lines 
A similar behaviour was long ago reported by Fetter JU| and Sonin ^2] for the low- 
frequency Kelvin modes, showing that the dissipation of such modes remains essentially 
equal to that of strictly rectilinear vortices. 

Our main objective in the present paper will be to analyze, within the framework of 
linear response theory the non-equilibrium dynamics and the equilibrium quantum 
fluctuations of the energy of the cyclotron modes. There exists an extensive literature 
on quantal Brownian motion of harmonic oscillators 14, , but it is important to realize 
that our problem presents a number of distinctive features that are not found in previous 
treatments, namely 

(i) Most of such previous studies have been focused on the coordinates or the 
momentum of the oscillator, rather than on the energy, eg treatments of the time 
correlation function of the energy are rather uncommon. 

(ii) It is evident that we are dealing with a very special heat bath, since in addition to 
being formed by Fermi particles and Bose quasiparticles, such bosons are characterized 
by a complex dispersion relationship which gives rise to different species (phonons, 
rotons R + and R~). 

(iii) The drag force on vortices that has been experimentally detected arise from 
scattering, thus we leave aside from our study the phonon radiation damping ITT|. 
We note that such a scattering interaction Hamiltonian is also unusual since it must 
be nonlinear in the heat bath operators. In fact, most of the Brownian motion models 
assume that the heat bath couples linearly to the harmonic oscillator, but in our case 
it is easy to realize that the scattering events must involve products of creation and 
annihilation operators of the particles that collide with the vortex. 

(iv) Our recent study [TT] has shown that the drag force could be affected by 
appreciable memory effects, so it will be important to extend our treatment beyond the 
usual Markovian approximation to explore such possible effects. 
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A suitable formalism to handle the above items can be found in reference [Tb\ , where 
a non-Markovian calculation of the energetic susceptibility of a harmonic oscillator, 
weakly coupled to boson and fermion environments, was carried out. So, we shall base 
our treatment on the above formalism. 

This paper is organized as follows, in the following section we describe our quantum 
model for the dissipative vortex dynamics which leads to an analysis of vortex oscillations 
and memory effects. Next in section 3 we summarize the main results of linear response 
theory applied to the vortex energy. In section 4, based on the previous calculation 
of the harmonic oscillator susceptibility, we analize memory corrections to the Markov 
approximation. In section 5 we calculate the response and time correlation functions, 
and study the dependence of the relaxation frequency on temperature and impurity 
concentration. Finally, in section 6 we gather the summary and main conclusions of our 
study. 



2. Quantum model for vortex dynamics 

We start from a vortex Hamiltonian given by, 

tf^) = ^(v 2 + ^-r 2 ), (9) 

where 

v = p/m„ + -zxr (10) 

corresponds to the vortex velocity r and p denotes the vortex canonical momentum. 
The second term in (110)1 corresponds to that of the vector potential (central gauge) in 
the electromagnetic analogy, and the z dependence in (JHJ) which arises from r(z) and 
p(z), corresponds to the rotation in the x-y plane parametrized by z, which results 
from following the helix path. Note that both the canonical momentum p and the 
Hamiltonian Q are given per unit length of the z axis. Then, it is easy to verify that 
the Hamilton equations lead from (jHJ) to the equation of motion (jSJ). 

The two-dimensional coordinate r of the vortex core can be written as the sum 
of the center coordinate r of the cyclotron circle plus the relative coordinate r' from 
such a center. Then, the quantization of such variables straightforwardly arises from 
the electromagnetic analogy: 

1 



r„ n r [(pt + /3)x + i(fl-P)y] (11) 

[(a f + a)x + i{a - a))y], (12) 



1 



where (/3') denotes a creation operator of right (left) circular quanta ^H] and L 
denotes the vortex line length. The z dependence of r arises from the replacements 
ft — ► exp(— ikz) fll, a) — ► exp(ikz) a\ in (JTT|l and (JUJ), and correspondingly for 
the annihilation operators. Actually, Fetter's theory identifies yq(z) as the whole 
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L 



displacement from the z axis (see ^U]j equation (14)). Analogous quantization for 
the canonical momentum p leads through (JTUJl and Q to a vortex Hamiltonian 

dzHJz) = h{n + uVa{a k + ±) + huj_{(3lf3 k + §) + ?koJatpt + a k p k ), (13) 

o 

where it is worthwhile noticing that v 2 and r 2 in Q turn out to be independent of z, 
as expected. The above Hamiltonian can be written to first order in u_/Q as, 

m(ala k + l)+huj-{ptp k + %) (14) 

where both polarizations (cyclotron and Kelvin modes) become decoupled. To prove 
this approximation, we first note that the set of eigenfunctions of (J 14)) are represented by 
wave functions corresponding to well-defined values of both numbers of circular quanta, 
right and left [TS|. On the other hand, the Schrodinger equation for the Hamiltonian 
(fT3~j) . can be easily solved by noting that the term proportional to r 2 in the expression 0, 
can be added to the corresponding term arising from v 2 , yielding a Schrodinger equation 
formally equivalent to the one with u;_ = 0, whose solution is well-known. Thus, we 
find that to the first order in oj_/VL we obtain the same spectrum of eigenvalues and 
eigenfunctions as from ()14)) . except for a slight correction in the radial coordinate of the 
wave functions, which has to be multiplied by the factor 1 + U-/Q. 

The Hamiltonian ()14)) corresponds to helical oscillations of fixed wavelength A = 
2n/\k\. In the final step of our quantization procedure we shall assume that the system 
obeys periodic boundary conditions over a length L along the z axis, so k will be 
restricted to values 2ns/ L, where s is a positive or negative integer. Thus, the complete 
vortex Hamiltonian is obtained by summing up the expression (|T4^ over all these values 
of k: 

H V = J2 hQ(k)(ala k + §) + hw-(k)(J3i/3 k + |). (15) 

k 

The above Hamiltonian differs from the one of Fetter's theory by the presence of the 
cyclotron modes of frequency Q(k) ~ Q (cf [TQJ , equation (11)). 
The heat bath Hamiltonian is given by, 

H B = Y^ b{. 6 k + £q C qj(T C qj(J , (16) 

k q,^ 

where 6 k denotes a creation operator of 4 He quasiparticle excitations of momentum 
fok and frequency u> k , and c qcr denotes a creation operator of solvated 3 He atoms of 
momentum frq, energy e q and spin 1/2 projection a. Note that we disregard any 
interaction between the heat bath particles themselves, since we shall work at low enough 
temperature and impurity concentration, so that such particles remain dilute allowing 
their treatment noninteractine; gas. 

To model the scattering interaction Hamiltonian, we will consider a generic 
momentum-conserving form: 

\ L dz £ [Ag bl 6 q + Tg 4 )CT c q , CT ]e-^)- r , (17) 

k,q,o- 
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where and r^ q denote scattering amplitudes depending on the momentum of the 
heat bath scatterers and the wave vector kz of the vortex wave. Recalling that the 
vortex coordinate can be written as r = r (z) + r'(z) + zz and taking into account 
that i"o(z) and r'(z) commute, the exponential factor in (fTTj) can be factorized as 
e -i(fc«-g»)« e -i(k-q)-r («) e -t(k-q).r (z)_ gj nce the amplitude of the vortex wave was assumed 

to be very small, it is tempting to expand the last two exponentials retaining only 
first order terms in r'(z) and r (z). This procedure was analyzed by Fetter ^U] for 
ro(2), finding that it leads to divergences at long wavelengths. The physical reason 
for this result can be understood by recalling that Vq(z) is linear in the creation and 
destruction operators, f3\ and f3 k . Consequently, an expansion in powers of r (z) is bound 
to fail whenever the energy per quantum huj- becomes very small at long wavelengths, 
as the transitions should involve many of these "soft" quanta ^U]- Notice that this 
argument does not apply to r'(z) since the operators a\ and at correspond to the high- 
frequency cyclotron quanta. This means that the treatment of Kelvin modes represented 
by 1*0(20 turns out to be considerably more complicated than that of the cyclotron modes 
represented by r'(z). In fact, only the phonon drag force arising from ro(z) could be 
analyzed by Fetter JU], but we shall see that all sources of friction acting on r'(z) can 
be studied. To this aim, let us set tq(z) = in (|T7j) while retaining only the first order 
term in r'(^). Then, using the second-quantized expression for r'(z) (cf (|12p). performing 
the integral in z and recalling the above-mentioned periodic boundary conditions, the 
interaction (fTTj) reads 

L E ( A kq ^ K + r kq Cq,«r] {<**„,, + 1^ jA K^X ~ k x ){h,k z -qA + S k ,q z -k z a k ) 

+ (Qy - k y)( 6 k,k z - q A - 5 kAz - kz a k )}}, (18) 

where the Kronecker-delta factors represent ^-momentum conservation and the 
scattering amplitudes were assumed to be independent of the sign of k, ie the interaction 
should be the same for both possible directions of a helical deformation. Note that the 
first term between braces in (J18)) does not contribute to the interaction, so it should be 
added to the heat bath Hamiltonian (fTtjj). Finally, summing up the expression (JTHJ) over 
k we obtain the interaction Hamiltonian: 



Him = Jt l - E [A k fc q b{ + r£-*> c q , CT ] 

V 47r Ps k q (7 

x[i(qx - k x )(a\ z _ qz + £%_ fc J + (q y - k y ){a\ z _ qz - a,,_*.)]. (19) 

From (|15|) and (|19|) . we may realize that each cyclotron mode of unperturbed 
Hamiltonian 

H k = hn(k)(ala k + l), (20) 
will evolve independently, interacting with the heat bath through the following terms of 

" (+) 



H itt = \j \ J2 t A kq h l K + r kq Ck,a Cq,«r] [fay ~ fcy) + *(<?* ~ k x )]a k 



■ 'kq "k "q r J kq ' k.rr ' q-"! V-'ln " « ■' ' " ' ' ' 'I • ".r/J" 

k,q,(T 
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(-) ) 
+ E [ A S 4 & q + rg Cq , CT ] [fo - q y ) + i(q x - k x )]a k , (21) 

k,q,cr J 

where the (±) sign above each summation symbol indicates that only the terms with 
k z — Qz — ±& must be considered. We recall, however, that for long wavelengths 
we have k <C a -1 ~ 1 A -1 , so it will be valid to neglect k in such z-momentum 
conservation relationships, except at extremely low temperatures (hc s k z ~ ksT). Thus, 
Hfo t becomes, 



' ^ s k,q,<r 

+ [(^ ~ + KQx ~ k x )]a k } (22) 

and the time evolution of a\ will be ruled by the Hamiltonian given by the sum of (J2(J|) . 
(fTp^) and ()22|) . To study the time evolution of the vortex coordinate r' = (x',y'), it will 
be convenient to use a complex form R' — x' + iy' since 

R'i^^Ze^al (23) 

is simply written as a linear combination of a[(t). In reference [TT] we studied the 
cyclotron dynamics of a rigid rectilinear vortex, this being equivalent to considering a 
single term with k —>■ in (|2*Hjl . We derived, within a weak-coupling approximation, a 
non-Markovian equation of motion for the mean value of the vortex position operator, 
finding that cyclotron frequency values within the range 0.01 — 0.03 ps -1 lead to a 
very good agreement with the experimental determinations of the longitudinal friction 
coefficient D (equation (j2J)), versus temperature and 3 He concentration. We showed that 
memory effects could represent up to ~ 10% of the D value as the number of heat bath 
scatterers is increased, that is, such effects are found to be increasing with temperature 
and impurity concentration. The scattering amplitudes leading to such results reads as 

HUE!, 



where V denotes the volume of the system, oj' k denotes the quasiparticle group velocity 
and <7o=18.54 A, corresponds to an effective cross section for vortex- 3 He scattering. Note 
that these amplitudes are in fact negligible with respect to the heat bath single-particle 
levels, LA^ ^ hu k and e k , for experimental sizes jTH] and not extremely low 

temperatures. 

The thermal excitation of vortex waves can be made consistent with the above 
experimental data, if we assume that the scattering amplitudes and in (|2*2|) 
are well approximated by the k = values, and (j2SJ), respectively. Note that 
this approximation is similar to the previous one, k z — q z = ±k — > (below equation 
f|2*T]) ) and also to Q(k) ~ Q in f|2T)j) . Thus, each a k (t) in (j2l?|) will present the same 
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dissipative evolution as a (t) ie, the same friction coefficient D should be ascribed to 
all long- wavelength cyclotron modes. This generalizes the previous result JOJ^], that 
the phonon friction coefficient associated to low-frequency Kelvin modes turns out to 
be essentially the same as that of strictly rectilinear vortices. 



3. The vortex energy in linear response theory 

According to the standard framework of linear response theory we will assume that, 
having the vortex reached thermal equilibrium with the heat bath before t = 0, a weak 
perturbing time dependent scalar field X(t) is coupled to the vortex Hamiltonian from 
t — onward. Then, the Hamiltonian of the whole system can be written, 

Hit) = H V + H B + H mt - X(t)H v , (26) 

where Hb is given by (|16j). the vortex Hamiltonian is given by, 

H v = mY t (fiUk + \) (27) 
k 

and the interaction Hamiltonian is given by, 



Hint = J 5 k*,q* t A kq & k + IfaJ 4,<T C <iA { [(?» ~ K) + KQx ~ ^)]4 

+[{k y - q y ) + i(q x - k x )]a k }. (28) 
Then, the mean value of the vortex energy can be written to the first order in X(t) as, 

(H v (t)) = (H v ) eq + f dr\(t - r)a(r), (29) 



where (H v ) eg = Nhfl{[exp(hfl/k B T) — + corresponds to the canonical equilibrium 
value (N = total number of long- wavelength cyclotron modes), and the function a(r) 
embodies the vortex response to the applied field. In particular, for a Dirac delta impulse 
X(t) = r 5(t — t ) , the above equation yields, 



(H v (t)) - (H v ) eq 



a{t - 1 ) (30) 



TO 

that is, a(r) represents the energy displacement from the equilibrium value, per unit 
strength of a pulse acting at r = 0. Now, if we assume a constant field X(t) = A , the 
so-called response function [TSJ is given by, 

= lim [(H v (t)) - (H v ) eq ]/X = fdra{r). (31) 

Ao^O Jo 

Finally, if the field is oscillatory X(t) = Aocos(a;£) (t > 0), the nontransient regime JH] 
can be described by setting t = oo in the upper limit of the integral in (|2U|) . 

roo 

(H v (t)) NT - (H v ) eq = X drcos[uj(t - r)] a(r) 

Jo 

= Aq Re [a (a;) exp(— iuit)] (32) 
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where a(u) may be defined as a complex generalized susceptibility, which is given by the 
Fourier-Laplace transform of the pulse response a(r), 



a(uj) = / exp(iur) a(r) dr. (33) 
J o 

Then, according to (J3~3~|) and (}3*T]) . the static susceptibility a(0) is given by, 



x Ninny 

a( w 0) = (i -> oo) 



sinh 



2fc B T 



(34) 



4fc B T 

where the right-hand side arises from taking into account that (H v (t — *> oo)) in (|31|) 
corresponds to the canonical distribution of a vortex with a Hamiltonian (1 — Xq)H v (or 
equivalently, a vortex with the Hamiltonian H v at the effective temperature T/(l — Ao)). 

As a final remark we note that from the susceptibility a (a;), one can readily get the 
equilibrium time correlation function, 

C(t) = ^(H v (t)H v (0) + H v (0)H v (t)) eq - {H v )% (35) 
via the fluctuation-dissipation theorem jT3j : 

C(u) = /icoth (^Lj lm[a(oj)], (36) 

where C(u) = /f^ dtexp(zu;t)C(t) denotes the Fourier transform of C(t). 



4. Analytic continuation of the generalized susceptibility: study of memory 
effects 

Being a Laplace transform, the generalized susceptibility (}33|) can be regarded as a 
function of a complex variable z, a(z), which must be analytic in the upper half-plane, 
Im z > 0. The important information, however, lies in the lower half-plane, where 
the spectrum of singularities of its analytic continuation yields the set of characteristic 
frequencies in the time evolution of the pulse response a(r) [T^j. Our calculation of 
a(z) is completely analogous to the one leading to the harmonic oscillator susceptibility 
in reference ^3]. Thus, we refer the reader to that paper for the technical details, and 
only quote here the final result (cf [T3], equation (2.25)): 

fit A = Nh 2 tf[q(z)-q(0)]/z 
K) [l-exp(-hn/k B T)][z + iv(z)Y { ' 

where q(z) and v(z) are Cauchy integrals, 

q(z) = / F(u) 38 

Hy > 2ni J-oouj - z y ' y ' 

viz) = / u(u), 39 

y ' 2ni J-oou - z y ' y ' 

with kernels (cf [T3], equations (2.31) and (2.32)), 

R(u)n(n + u) R(-u)n(n - u) 

t [LO) = — r (4UJ 

n(uj) n{—uj) 
u(u) = -[R(u) + R(-lu)], (41) 

% 
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being, 



n{ui) = [exp(hu /ksT) — 1] 



-i 



(42) 



and 



R(u) 



2% 



(43) 



The function D in the above equation has been studied in previous works ITT] , 
since D(Q) corresponds to the longitudinal friction coefficient in the Markovian 
approximation. It reads jTTj, 



where n k = [exp(hu k /k B T) — 1] _1 and f k = {exp[(efc— /x)/A; B T] + l} _1 respectively denote 
the thermal equilibrium Bose and Fermi occupation numbers for the corresponding 
scatterers. The expression (pT7jl for the susceptibility, on the other hand, is fully non- 
Markovian and each pole Zj of it yields a term proportional to exp(— %Zjj) in air). 
Then, the Markov approximation consists in neglecting the set of such poles which are 
located far enough from the origin, so that they yield rapidly vanishing terms, ie terms 
which decay faster than any observational timescale. This is the case for the set of 
poles arising from q(z) in (J37)l . actually poles of n(fl ± uS) in (|3Uj) [To] , which are of 
the form ±Q — in2nkBT/h (n =1,2,...). Such poles give rise to exponentially decaying 
terms in the expression of a(r), which have lifetimes shorter than h/ksT. This thermal 
timescale turns out to be much smaller than the one arising from the friction coefficient, 
ie h/k B T < m v /D(Q) = [p^/DiQ^Q- 1 for T < 1.5 K and fi's of order 0.01 ps" 1 . So, 
the above set of thermal poles can be safely ignored. It is clear then, that we should 
look for poles of order Q[D(Q)/ p s h] and, in the following, we shall see that they arise 
from the equation z + iv(z) = 0. In fact, looking for a solution close to the origin, one 
may begin with the ansatz z = —iv(0) and next proceed iteratively, ie Z\ = —iv(z ) 
and so on. This leads to a very rapid convergence to a solution z s that is better worked 
out in terms of the Taylor expansion of v(z) around the origin: 



where the second and third term inside the square brackets represent first and second 
order corrections to the zeroth order solution, respectively. The Cauchy integral 
and its derivatives in ([45)1 . can be written as [15], 



D(Q) 




(44) 



Zs = -iv(0)[l - iv'(0) - v"(0)v(0)/2] 



(45) 



v(0) 



f(0)/2 = 2VLD{VL)/p s h 



(46) 




i7rp s h Jo u 
2QD{Q)} 



(47) 
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v "(0) = z/"(0)/2 = ^-YlD'iSl) + VLD"(n)\. (48) 
p s h 

We see from fl4*U|) that the zeroth-order solution z = —iv(0) has in fact the expected 
dependence and, moreover, the factor 2 in the expression of v(0) is easily interpreted if we 
recall that a damping in the velocity like exp(—QtD(Q)/p s h) should give rise to a twice 
faster energy damping. The first and second order corrections in (|43j) are easily evaluated 
from (J47|) and (J48|) . and they are always negligible, eg for ordinary helium at T = 0.67 
K (Q = 0.01 ps" 1 ) we have -iv'(0) = 1.52 x 10~ 5 and -v"(0)v(0)/2 = 1.32 x 10" 9 . 
In conclusion, we have found that the Markovian pole —i2QD(Q)/p s h is unaffected by 
memory corrections. This is to be contrasted with the appreciable memory corrections 
to the friction coefficient seen in section 2. Therefore, the Markovian approximation for 
the generalized susceptibility (|3"Tj) reads, 

Nh 2 Q 2 q'{0) 

aM[z) ~ [1 -ex V (-hQ/k B T)][z + iv(™ (l9) 



which should be valid for z inside a circle with a radius r fulfilling v (0) < r k B T /h. 

Notice that we have replaced the expression [q(z) — q(0)]/z in (j37J) by the derivative 
q'{0). This approximation may be readily tested if one considers the first order term 
q"(0)z/2 at z = z Q = — iv(0). In fact, we have [T3] 

g '(0) = F'(0)/2 = -n(n)z Q /k B T (50) 

and 

</'(0) = - r-AF{u)-F'{Q)ul (51) 



m Jo uj 3 

so, we may find again that the first order correction is totally negligible, eg for ordinary 
helium at T = 0.67 K (Q = 0.01 ps -1 ) it represents ~ 10 -6 of the zeroth order ^'(O). 
Then, replacing ()50)1 in |i9|) we have the final expression, 

Nh 2 n 2 i2nD(n)/ Ps h 

4k B T smh 2 ( y hQ/2k B T) [z + i2QD {tt)/p s h] 
which immediately reproduces the result (|3^j) for the static susceptibility a(0). 



a M{z) = - ^ 2^r^ /»,. ^ r. , t~\ / r^\\ / , ,,n ( 52 ) 



5. Response and time correlation functions 

From (|52j) one easily extracts ol(t) and the response function (|31l) . 



-r / \ N(my 



4k B T 



sinh 



2k B T / 



1 -2 



[1 - exp(-2tttD(n)/p a h)]. (53) 



This kind of response, characterized by a single relaxation time, is well-known in theories 
of dielectric and magnetic relaxation and goes under the name of Debye response |13j . 
The Fourier transform of the time correlation function arises from (|36|) and (J52j) : 



ik B T 



sinh 



2k B T i 



Ticoth.(w)we , , 

2,2 ' 54 
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where e = [QD(Q)/ p s h}/[k B T/h] and w = hu/2kBT. Being e < 1, the function of w, 
ej (e 2 + w 2 ) in (j54j) turns out to be sharply peaked around w = 0, so we may approximate 
coth(u>)u> ~ 1 and thus get the antitransform: 



Therefore, the time correlation function is ruled by the same relaxation time of the 
response function, as expected. Notice that C(0) in (|B3j) actually corresponds to 
(H 2 ) eq — (H v )l q , as can be easily verified by an elementary calculation of (H 2 ) eq in 
the canonical ensemble. 

To conclude, it is interesting to analyze how the relaxation frequency 2QD(Q)/ p s h 
depends on temperature and impurity concentration. In figure 2, D(Q)/p s h is plotted 
against temperature for several 3 He concentrations [3]. Such curves actually correspond 
to O = 0.01 ps -1 , but it is important to remark that the dependence on Q turns out 
to be negligible for JTs within 10~ 2 ps -1 or less JZIE]- It is also worth noticing that 
consistently with our weak-coupling approximation jT3 El CHj , the Markovian friction 
coefficient D(Q) always remains small compared to the coefficient p s h of the Magnus 
force in (JTJ) (actually p s h has virtually no dependence on temperature for T < 1 K). 
The lowest curve in figure 2 corresponds to pure 4 He and it displays two well separated 
regimes [T7| . 



where A/fcg=8.62 K corresponds to the height of the roton minimum in the dispersion 
curve of figure 1. Note that the phono n-roton transition clearly manifests itself as an 
intermediate region of positive second derivative. The remaining curves in figure 2 
correspond to finite 3 He concentrations that have been experimentally studied [3] • Such 
concentrations are low enough to allow a Maxwell-Boltzmann approximation for the 3 He 
statistics in (f4~4*j) . Then, the low-temperature regime of D(Q), which is now dominated 
by impurity scattering, turns out to be proportional to y/T and 3 He concentration [*Hlllj. 
We may see from figure 2 that phonon effects are completely hidden in ordinary helium, 
since the 3 He domain extends as far as T ~ 0.4 K. For higher concentrations such 
domain reaches higher temperatures hiding also the first portion of the roton curve. 

Finally, it is worthwhile observing that relaxation frequencies with any temperature 
dependence, rarely appear in models of quantal Brownian motion of harmonic oscillators, 
since most of them assume, in contrast to our scattering model, linear couplings in the 
heat bath operators [14*1 115]. 

6. Summary and conclusions 

A generalization of the quantum theory of vortex waves [""J [TO], has been proposed 
to study the damping of cyclotron modes. We have shown that the friction values 
should be practically unaffected by such oscillations, a result which was already known 



C(t) 



N (hn) 2 exp[-2fl|t|D(fi) I p a h] 
4sinh 2 ("in/2fc B T) 



(55) 




T 5 , for T < 0.4 K (phonon domain) 

exp(— A/ksT), for T > 0.5 K (roton domain), 




Figure 2. Relative value of the friction coefficient D(£l) with respect to the Magnus 
force coefficient p s h, versus temperature for several 3 He concentrations C. 

in the case of phonon scattering of low- frequency Kelvin modes [TTH H2] . All sources 
of dissipation arising in ordinary helium, viz phonons, rotons and 3 He atoms were 
considered, showing that appreciable memory effects must be taken into account in the 
evaluation of the friction coefficient. We have also analyzed memory corrections to the 
Markov approximation in the case of the nonequilibrium energetics of cyclotron modes, 
finding this time that they are negligible. We have shown that the vortex response 
is of the Debye type, ie it is ruled by a single relaxation frequency which governs the 
time correlation function as well. Such a relaxation frequency is shown to embody all 
the complexity of the heat bath, in that very well separated regimes belonging to the 
different species comprising the superfluid helium, are recognized from its dependence 
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on temperature and impurity concentration. 
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